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We search for B° meson decays into two-body combinations of K° , n, r/' , and cj> mesons in 324 
million BB pairs collected with the BABAR detector at the PEP-II asymmetric-energy e + e~ collider 
at SLAC. We measure the following branching fractions (upper limits at 90% confidence level) 
in units of 1CT 6 : B(B° -> nK a ) = 1.81^ ± 0.1(< 2.9), B(B° -> ryq) = l.ltg;* ± 0.1(< 1.8), 
B(B° -> ri(t>) = 0.1 ± 0.2 ± 0.1(< 0.6), B(B° -> ij» = 0.2±g;| ± 0.1(< 1.0), and B(B° -> n'r,') = 
l.Olo.e ± 0- 1 (< 2 - 4 )> where the first 

error is statistical and the second systematic. 
PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



We report the results of searches for B° or B° me- 
son decays to two charmless pseudoscalar mesons Q 
r]K°, rfrf , and to the pseudoscalar- vector combina- 
tions r]4>, r/'(f). None of these decays has been observed 
previously; the published experimental upper limits on 
their branching fractions lie in the range (2 — 10) x 1CP 6 
0, The theoretical predictions for these branching 
fractions are less than a few per million by most estimates 
H H H HH H E3 Theoretical approaches include those 
based on flavor SU(3) relations 0,113, EJ, effective Hamil- 
tonians with factorization and specific B-to-light-meson 
form factors Q, perturbative QCD QCD factoriza- 
tion 0, and soft collinear effective theory (SCET) |Toj | . 
Important advances in the theoretical understanding of 
hadronic charmless two-body B meson decays have oc- 
curred in the past few years With more precise ex- 
perimental results one can test and constrain the models. 
Improved measurements of decays with isoscalar mesons 
can also help to better understand the large difference be- 
tween the branching fractions for B — ► rf ' K and B — > r\K 
decays [HE3- 

Branching fractions or limits in the 7777, rfrf, 770, and 
7/0 channels are relevant for the accuracy with which CP- 
violating asymmetry measurements can be interpreted. 
The coefficient S of the CP- violating sinusoidal factor in 
the time evolution of ij'K and 4>K° can be related to 
the CKM phase (3 = arg {-V cd V* b /V td V? b ) if these de- 
cays are dominated by a single weak phase 0]. Addi- 
tional higher-order amplitudes with different weak phases 
would lead to deviations AS* between the value measured 
in these rare modes and the precise determination in the 
more copious B° decays to charmonium-if final states. 
SU(3) flavor symmetry [Til fl5j | relates the strength of 
such additional amplitudes to the decay rates of certain 
two-body B decays, including 7777, T}'r]' , rj(j), and rj'<j>. 

The results presented here are based on data collected 
with the BABAR detector at the PEP-II asymmetric- 
energy e + e~ collider located at the Stanford Linear Ac- 
celerator Center. An integrated luminosity of 289 fb _1 , 
corresponding to N BB - — 324 million BB pairs, was 
recorded at the T(4S) resonance (center-of-mass energy 
y/s= 10.58 GeV). 

Charged particles produced in e + e~ interactions are 
detected, and their momenta measured, by a combination 
of a vertex tracker, consisting of five layers of double- 



sided silicon microstrip detectors, and a 40-layer cen- 
tral drift chamber, both operating in the 1.5 T magnetic 
field of a superconducting solenoid. We identify photons 
and electrons using a CsI(Tl) electromagnetic calorime- 
ter. Further charged-particle identification is provided by 
the average energy loss (dE/dx) in the tracking devices 
and by an internally reflecting ring-imaging Cherenkov 
detector (DIRC) covering the central region. 

We select 77, 77', 4>, p°, K s , and 7r° candidates through 
the decays 77 -> 77 (?7 77 ), 77 -» 7r+7r~7r° (773^), 77' -> 
7771+71-- with 77 -> 77 (77^), 7/ -> p°7 {rfpj, 4> -> 
K+K-, p° -> tt+tt-, K s -> tt+tt-, and 7T° 77. The 
photon energy must be greater than 30 (100) MeV for 
7T (prompt 77 from B) candidates, greater than 200 MeV 
in 77' — * fry, and greater than 50 (100) MeV in 77' w (in 
the B — * 77' w r]' V7T7T decay mode). We make the follow- 
ing requirements on the invariant masses (in MeV/c 2 ): 
490 < m 77 < 600 for 77 77 , 120 < 777 77 < 150 for tt , 
510 < m„ n < 1000 for p°, 520 < m vww < 570 for 773^, 
930 < 777^^^ < 990 for 7^ )7r7r , 910 < m PJ < 1000 for 77^, 
1005 < m, K+K - < 1035 for </>, and 486 < 777^ < 510 
for Kg. For K s candidates we also require a vertex 
X 2 probability larger than 0.001 and a reconstructed de- 
cay length greater than three times its uncertainty. Sec- 
ondary charged pions in 77 and 77' candidates are rejected, 
if their DIRC and dE/dx signatures are consistent with 
protons, electrons, or kaons. Similarly, tracks from <f> 
decays are required to be inconsistent with protons, elec- 
trons, and pions. 

A B meson candidate is characterized kinc- 
matically by the energy-substituted mass toes = 
[(^s + p • Pb) 2 /Eq — p 2 B ]? and energy difference 
AE = E* B — \\fs, where the subscripts and B refer to 
the initial T(4S) and to the B candidate, respectively, 
and the asterisk denotes the T(4S) rest frame. 

Backgrounds arise primarily from random combina- 
tions of tracks and neutral clusters in e + e~ — > qq contin- 
uum events, where q = u,d, s, or c. We reject these events 
by using the angle £?t between the thrust axis of the B 
candidate in the T(AS) frame and that of the rest of the 
event. The thrust axis of the B candidate is obtained as 
the thrust axis of the B decay products. The distribution 
of |cos6*t| is sharply peaked near 1.0 for combinations 
drawn from jet-like qq pairs, and is nearly uniform for 
T(45) -> BB events. We require | cos6» T | < 0.9. To dis- 
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criminate against r-pair and two-photon backgrounds we 
require the event to contain at least three tracks or one 
track more than the topology of our final state, whichever 
is larger. In decays containing a prompt f? 77 from B we 
require \H n \ < 0.9 to remove random combinations with 
soft photons, where H v is defined below. If an event has 
multiple B candidates, we select the candidate with the 
highest B vertex \ 2 probability or using a \ 2 quantity 
computed with the rj or r\' masses, depending on the de- 
cay mode. More details on the analysis technique can be 
found in Ref. [l7| . 

We obtain yields from unbinned extended maximum- 
likelihood (ML) fits. The principal input observables are 
AE, toes, and a Fisher discriminant T Where rel- 
evant, the invariant masses m ros of the intermediate res- 
onances and angular variables Ti. defined below are used. 
The Fisher discriminant T combines four variables: the 
angles with respect to the beam axis of the B momen- 
tum and B thrust axis (in the T(45) frame), and the 
zcroth and second angular moments io,2 of the energy 
flow about the B° thrust axis. The moments are defined 
by Lj = J2iPi x \cos 9i\ J , where 9% is the angle with re- 
spect to the B thrust axis of track or neutral cluster i, 
Pi is its momentum, and the sum excludes the B candi- 
date. For ry 77 (0), Ti. v (He/,) is defined as the cosine of 
the angle between the direction of a daughter 7 (K) and 
the flight direction of the parent of 77 (0) in the rj (</>) 
rest frame; for rj' , 7i p is the cosine of the angle between 
the direction of a p daughter and the flight direction of 
the rj' in the p rest frame. The set of probability density 
functions (PDF) used in ML fits, specific to each decay 
mode, is determined on the basis of studies with Monte 
Carlo (MC) simulated samples [l^. We estimate BB 
backgrounds using MC samples of B decays. The esti- 
mated BB background is found to be negligible for all of 
our decay modes except rj^^Kg and ?? 77 0. 

The extended likelihood function is 



N 



£ = exp(-^rij) Y[ 



i=l 



3 = 1 



(1) 



where N is the number of input events, rij is the number 
of events for hypothesis j (J — 1 for signal, j = 2 for 
continuum background, and j = 3 for BB background), 
and Vj(x.i) is the corresponding PDF evaluated with the 
observables Xj of the i th event. The BB background 
component is used in the decay modes rj^Kg and ?7 77 </>. 
Since the correlations among the observables in the data 
are small, we take each Vj as the product of the PDFs for 
the separate variables. We determine the PDF param- 
eters from simulation for the signal and from sideband 
data (5.25 < m ES < 5.27 GeV/c 2 ; 0.1 < \AE\ < 0.2 
GeV) for continuum background. We float some of the 
continuum PDF parameters in the ML fit. We parame- 
terize each of the functions 'Pi(toes): Vi(AE), Vj(F), 
and the peaking components of Vj(m les ) with either a 



Gaussian, the sum of two Gaussians, or a Crystal Ball 
function [2(J as required to describe the distribution. 
Slowly varying distributions (m rcs and AE for combina- 
torial background, and angular variables) are represented 
by linear or quadratic functions. The combinatorial back- 
ground in toes is described by the ARGUS function [2l| . 
Large data control samples of B decays to charmed final 
states of similar topology are used to verify the simulated 
resolutions in toes and AE. Where the control samples 
reveal differences between data and MC in mass or en- 
ergy resolution, we shift or scale the resolution used in 
the likelihood fits. The bias in the fit is determined from 
a large set of simulated experiments, each one with the 
same number of qq and signal events as in data. 

Table [H shows the measured yields, efficiencies, and 
products of daughter branching fractions for each decay 
mode. The efficiency is calculated as the ratio of the 
numbers of signal MC events after the cut based selection 
to the total generated. We compute the branching frac- 
tions from the fitted signal event yields, reconstruction 
efficiency, daughter branching fractions, and the number 
of produced B mesons, assuming equal production rates 
of charged and neutral B pairs at T(AS). We correct the 
yield for any bias measured with the simulations. We 
combine results from different channels by adding the 
values of — 21n£ (parameterized in terms of the branch- 
ing fraction), taking into account the correlated and un- 
corrected systematic errors. We report the statistical 
significance and the branching fractions for the individ- 
ual decay channels. For the combined measurements we 
also report the 90% confidence level (CL) upper limits. 

The statistical error on the signal yield is taken as the 
change in the central value when the quantity — 21n£ 
increases by one unit from its minimum value. The sig- 
nificance is taken as the square root of the difference be- 
tween the value of —2 In £ (with systematic uncertainties 
included) for zero signal and the value at its minimum. 
We determine a Bayesian 90% CL upper limit assuming 
a uniform prior probability distribution by finding the 
branching fraction below which lies 90% of the total of 
the likelihood integral in the positive branching fraction 
region. 

Figure^shows, for representative fits, the projections 
onto Toes and AE for the five decay modes. The points 
show the data after a channel-dependent requirement on 
the probability ratio V\/(V\ + Vn + V3), optimized to 
enhance the signal sensitivity and with the probabilities 
Vj evaluated without using the variable plotted. The 
solid curves show the total rescaled fit functions. 

The main sources of systematic error include uncer- 
tainties in the PDF parameterization (0-2 events) and 
ML fit bias (0-2 events) . We evaluate these uncertainties 
with simulated experiments by varying the PDF param- 
eters within their errors and by embedding MC signal 
events inside background distributions simulated from 
PDFs. The uncertainty on N B g is 1.1%. Published world 
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TABLE I: Fitted signal event yield, fit bias, detection efficiency e, daughter branching fraction product significance S, 

and measured branching fraction B with statistical error for each decay mode. For the combined measurements we give the 
significance (with systematic uncertainties included) and the branching fraction with statistical and systematic uncertainty (in 
parentheses the 90% CL upper limit). 



Mode 


Yield (ev) 


Fit bias (ev) 


e (%) 


WBi (%) 


S (a) 


/3(10 -6 ) 






1Q + 10 

t»_9 


i n o in c 
+U.o ± U.D 


zd. / ± U.y 


1Q C 
lO.D 


Z.U 


1 c+0.9 
1 -°-0.8 






1 1+6 


i i i i n a 
+ 1.1 ± U.4 


1 f .O ± U.D 


~7 Q 
( .O 


O *7 
Z. / 


2 - 4 -l.l 




TjK 










3.5 


1.8I 6 ± 0.1 


(< 2 - 9 ) 


V-rt V-n 


17 -9 


+o.y ± U.D 


on piio 
zU.o ± l.o 


lO.D 


i.y 


1 o + l.O 
± -°-0.9 






m+ 7 
10 -5 


+U.D ± U.4 


1 O O _l_ 1 o 


i t n 

i f .y 


O 1 

2.1 


n n+0.6 
0-9-0.5 






o+3 
2 -2 


+U.O ± U.4 


11. D =t U.o 


. 1 


i i 
1.1 


1 1+ 1 ' 6 
i - 1 -1.0 




vv 










3.0 


1.1±°-J±0.1 


« 1-8) 


Vty<P 


-u±S 


-2.4 ±0.6 


32.3 ± 1.2 


19.4 


0.0 


-0.4±^ 




V3w<j> 




+0.8 ±0.3 


20.7 ± 1.0 


11.1 


1.5 


o.7±g:S 




ricj) 










0.0 


0.1 ± 0.2 ± 0.1 


« 0.6) 


Vp-y<t> 


i+3 

A -2 


-0.6 ±0.3 


23.1 ± 1.1 


8.6 


0.8 


u - J -0.3 




o+9 


-1.0 ± 0.4 


22.5 ± 0.9 


14.5 


0.0 


n n+0.9 














0.5 


0-2^ ±0.1 


« 1-0) 


T]ri TV TV Vl] TV TV 


l + -t 


±0.3 ±0.2 


15.2 ± 1.0 


3.1 


1.2 


0.8+J;? 




V^TTTrVp-'f 




±1.5 ±0.3 


17.6 ±0.8 


10.3 


1.5 


1 2+ 1 - 1 




rj'rj' 










1.8 


i.o±°:5 ±o-i 


« 2.4) 




m ES (GeV/c*) AE (GeV) 



FIG. 1: Signal enhanced projections on tties (left) and AE 
(right) in the decays: (a, b) 77^°, (c, d) 7777, (e, f) 77^, (g, h) 
rj'cj), (i, j) 7/77'. Points with error bars (statistical only) repre- 
sent the data (combined measurements), the solid line the full 
fit function, and the dashed line its background component. 

averages provide the uncertainties in the B-daughter 
branching fractions (1-7%). Other sources of systematic 
uncertainty are track (1-3%) and neutral cluster (2-6%) 
reconstruction efficiencies. The validity of the fit pro- 
cedure and PDF parameterization, including the effects 
of unmodeled correlations among observables, is checked 
with simulated experiments. 

Grossman el al. I1J| introduced a method to deter- 



mine a bound on |A5/| = \Sf — sin2/3| where / is 
a CP eigenstate produced in charmless B° decays and 
S is the coefficient of the CP-violating sinusoidal fac- 
tor mentioned above. The method relics on SU(3) fla- 
vor symmetry and the measured branching fractions of 
charmless, strangeness-conserving B° decays to constrain 
the unknown contributions of suppressed amplitudes in 
B° — > /. Two of the channels in our study, r\r\ and 
7/77', are relevant to the ASf bound for / = rj'K , while 
two others, r/cj) and 77' 0, are relevant for / = <f>K a . Us- 
ing the technique described in Ref. and evaluating 
90% CL upper limits, we find \AS V ' K o\ < 0.15 and 
\AS,f, K o\ < 0.38. This new AS^/^o bound also makes 
use of our recent results on the B° — > rj'rj, r/'n , and 
rjTT channels. 

In summary, we present updated measurements of 
branching fractions for five B° decays to charmless meson 
pairs. Our results represent substantial improvements on 
the previous upper limits [HQ 
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